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Spin-orbit interaction (SOI) leads to spin precession about a momentum- 
dependent spin-orbit field. In a diffusive two-dimensional (2D) electron gas, the 
spin orientation at a given spatial position depends on which trajectory the elec¬ 
tron travels to that position. In the transition to a ID system with increasing 
lateral confinement, the spin orientation becomes more and more independent on 
the trajectory. It is predicted that a long-lived helical spin mode emerges [T| 2]. 
Here we visualize this transition experimentally in a GaAs quantum-well struc¬ 
ture with isotropic SOI. Spatially resolved measurements show the formation of 
a helical mode already for non-quantized and non-ballistic channels. We find a 
spin-lifetime enhancement that is in excellent agreement with theoretical predic¬ 
tions. Lateral confinement of a 2D electron gas provides an easy-to-implement 
technique for achieving high spin lifetimes in the presence of strong SOI for a 
wide range of material systems. 

In a diffusive electron system with intrinsic SOI (e.g., of Rashba or Dresselhaus type), 
the effective spin-orbit field changes after each scattering event. This leads to a random¬ 
ization of spin polarization that is described by the Dyakonov-Perel (DP) spin-dephasing 
mechanism [3], in the case of an initially homogenous spin excitation. Given a local spin 
excitation, a spin mode emerges that is described by the Green’s function of the spin dif¬ 
fusion equation gHSj. For a 2D system in the weak SOI limit, analytical solutions exist for 
a few special situations, such as for the persistent spin helix case with equal Rashba and 
Dresselhaus SOI PHI]. In the isotropic limit (either only Rashba or only linear Dresselhaus 
SOI), the spin mode is described by a Bessel-type oscillation in space (see Fig. 3>)0- The 
spin lifetime of such a mode is only slightly enhanced [4] compared with the DP time because 
rotations about varying precession axes (see Figs, [lj^-e) do not commute and therefore the 
spin polarization at a given position depends on the trajectory on which the electron reaches 
that position. If the electron motion is laterally confined by a channel structure of width 
w, the spin motion is restricted to a ring on the Bloch sphere (see Figs. IH and ID 1 )- In this 
situation, the spins collectively precess along the channel direction (Fig. §) mm- This 
extends even into the 2D diffusive regime as long as the cumulative spin rotations attributed 
to the lateral motion are small, i.e., as long as wq° < 1, where q° is the lateral wave num¬ 
ber of the 2D spin mode. As a consequence, for a 2D diffusive system, increasing lateral 
confinement is predicted to result in an enhanced spin lifetime proportional to ( q°w ) 2 Q] 2 . 
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This effect could be highly relevant for spintronics applications because it circumvents the 
conventional trade-off between a long spin lifetime and strong SOI. It has been experimen¬ 
tally explored in different ways, including measurements of weak-antilocalization mm, 
the inverse spin-Hall effect [15], and time-resolved Kerr rotation |16| . None of these works 
were able to resolve the spin dynamics both spatially and temporally, and a quantitative 
investigation of the spin mode in the confined channel is still lacking. 

We experimentally explore the dynamics and spatial evolution of electron spins in a 2D 
electron gas hosted in a symmetrically confined, 12-nm-wide GaAs/AlGaAs quantum well 
where the linear Dresselhaus SOI is much larger than the Rashba or the cubic Dresselhaus 
SOI, thus providing an almost isotropic SOI. To study the transition from 2D to ID, we have 
lithographically defined wire structures along the [110] (x) and [110] (y) directions with the 
channel width w ranging from 0.7 to 79 /im. 

Figure [pi shows a sketch explaining the measurement principle. Spins polarized along the 
out-of-plane direction, z, are locally excited at time t = 0 by a focused, circularly polarized 
pump laser pulse, which has a Gaussian intensity profile of a sigma-width of 1.1 /im. A 
second, linearly polarized probe pulse measures the out-of-plane component, S z , of the local 
spin density using the magneto-optical Kerr effect. The spatial evolution of the spin packet 
is mapped out along the channel direction for various time delays, t , between pump and 
probe pulse. All measurements have been performed at a sample temperature of 20 K. 

A measurement of spatially resolved spin dynamics in a channel in the 2D limit (w = 
19 n m) is shown in Fig. [2^t. Spins are excited at t = 0 and at x = y = 0 and traced as a 
function of y and t. At y = 0, S z simply decays in time. It reverses its sign after t > 400 ps 
for electrons that diffused along y by more than 4 /an, seen as a faint blue color in the 
Figure. The situation is different in the 0.7-//m-wide channel (Fig. i>). Here, spin decay is 
strongly suppressed and S z reverses its sign multiple times along y at later times. Note that 
the pattern is overlaid with the spin texture that survived from the previous pump pulse at 
t = —12.6 ns. Figure^ shows measured data of S z (y) for the 19-/an and the 0.7-/an-wide 
channels taken at t = 1.5 ns. The comparison of the two curves clearly shows an enhanced 
S z and strong oscillations along y in the narrow channel. This indicates a helical spin mode 
in the ID case. The helical nature is further supported by measured maps where an external 
magnetic field is applied along the x direction, rotating the helix as a function of time, see 
supplementary information. 
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FIG. 1. Measurement principle and expected spin modes, a, A focused pump laser pulse 
locally excites out-of-plane spin polarization S z with a Gaussian width of ^1.1/im. A probe laser 
pulse measures the local spin distribution after a time delay t via the magneto-optical Kerr effect. 
Pump and probe spots are spatially scanned against each other, b, In the 2D case, spin diffusion in 
the isotropic SOI-field (e) turns a local spin excitation along z into a Bessel-type spin mode, whose 
S z component is shown (red color indicates positive, blue color negative S z ). The spin-trajectories 
in the 2D plane, c, are correlated with the trajectories on the Bloch sphere, d, such that the local 
spin density decays more slowly than that of the whole ensemble, f, in the ID case, the emerging 
spin mode is long-lived and described by a cosine oscillation of S z , corresponding to a helical 
rotation of the spin polarization, g, The lateral confinement restricts the diffusive trajectories in 
real space, such that the spins on the Bloch sphere evolve on a ring (h) whose width scales with 
the channel width, w. The smaller ic, the more the spins precess about a single axis, leading to a 
drastic increase of the lifetime of the helical mode. 

For a deeper analysis, it is advantageous to Fourier-transform S z (x,y,t) into momentum 
space. Thereby one obtains Fourier components S z (q X) q y , t) at wave numbers q x and q y that 
according to theory decay biexponentially in time [5j. For channels narrower than 15 /in 1 , the 
spin modes exhibit a pronounced structure only along the channel direction, and we therefore 
analyze the ID Fourier transformation along this direction. For wider channels, we obtain 
the 2D Fourier transformation from ID scans of S z by assuming a radially symmetric spin 
mode, see supplementary information for details. This is justified because we observe a 
similar dependence of S z along the x and y directions, as seen from the values obtained for 


4 

















FIG. 2. Direct mapping of spin precession and spin decay, a-b, Spatial maps along the 
channel direction y of the out-of-plane spin density S z for varying time delays t between the pump 
and probe laser pulses. The 19-^m-wide channel (a) represents a 2D situation, whereas the 0.7- 
/im-wide channel (b) is close to the ID limit and exhibits a long-lived mode with oscillating S z (y ) 
[sign encoded as red (+) and blue (-)]. The mode from the preceding laser pulse (pulse period 
12.6 ns) is still visible at negative delay, c, Line-cuts through the data of (a) and (b) at t — 1.5 ns 
for comparison, d and e show the Fourier transform of S z (y ) for the 19- and 0.7-/im channel, 
respectively. For w = 0.7 /im, the initially Gaussian spectrum quickly converges to a long-lived 
spin mode at finite q y = q y . f, Line-cuts through the data of (d) and (e) at q y ~ q y - At each q y , 
the amplitude decays biexponentially, with the the longer-lived mode visible at longer t. 


wavenumbers q® and q® later in the text. 

Figures[2]i and[2^ show S z (q y) t ) for the 19- and 0.7—gm wires, respectively. The Gaussian 
distribution of S z (q y) t) decays in time with very different rates for varying q y , which are 
minimal at a finite wavenumber, q Figure^ shows traces at q y ^ q y for the two cases. 
For t > 500 ps, we fit each trace with a single exponential decay to obtain the momentum- 
dependent lifetime r(q y ) of the longer-lived spin mode |5, [9j. The decay rates, l/r(q y ) are 
shown in Fig. |3^. In both the ID and the 2D case, 1/r vs q y can be well approximated close 
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FIG. 3. Fit results, a, Decay rates for the 19-ym (diamonds) and 0.7-/im-wide (dots) channels 
along the y direction as a function of the wave number, q y . Data is obtained from a fit of the 
Fourier-transformed S z (y) to an exponential decay for t > 500ps, thus corresponding to the long- 
lived mode. Solid lines are parabolic fits of the decay rate versus g, from which the position, q®, and 
the lifetime, r°, of the evolving spin mode are obtained, b, Values for the spin-mode wave numbers 
q® and q®, shown for measurements at various channel widths, re, along the x and y direction, c, 
Lifetime of spin modes, r°, as a function of w and for both channel directions. Solid black lines 
are the theoretically expected lifetimes. The yellow solid line is their interpolation. In b and c, 
circles (diamonds) stand for fit values obtained from ID (2D) Fourier transformations. 

to q^ by the parabolic function PE! 

1/t=1/t° + D s (q y -q° y ) 2 , (1) 

where D s is the spin diffusion constant [CJ. Figures [3p and §■ plot the values obtained 
for q [ y iq l l) and r°, respectively, for channels along the y ( x ) direction and of various widths. 

Comparing and q y in Fig. [iijr we observe a slight anisotropy characterized by > q y . 
This means that the SOI is stronger for electrons that move along x and indicates a remaining 
Rashba field due to a slight asymmetry in the quantum well. The SOI coefficients, a and /3, 
are obtained from q y and measured in the ID limit by using the expressions 


2m* 


h z 

qI = 


2 m* 

~ + P) 


0.7 ym 1 , and 

0.8 ym _1 . 
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( 2 ) 

( 3 ) 















Here, m* is the effective electron mass and h is the reduced Planck’s constant, a is the 
SOI parameter of the Rashba field and (3 = f3\ — /? 3 that of the Dresselhaus field. 3\ and /? 3 
characterize the linear and cubic Dresselhaus fields, respectively. Values for a, (3i, (3% and 
D s are given in the supplementary information. 

The dependence of q ]“ on w is rather flat, whereas decreases for increasing w. This is in 
agreement with the prediction that of the 2D spin mode is smaller for slightly anisotropic 
SOI than expected from Eq. © [51 H8]. Close to the persistent spin helix situation, the 
same effect leads to a suppression of precession along y. 

The lifetime, r°, however, behaves almost identically for both wire directions and increases 
by about one order of magnitude from w = 19 to 0.7 //in. Theory provides expressions for 
the lifetime in the 2D limit EEH3, r 2 d, and in the intermediate regime png i Tim, see 
supplementary information. For very narrow channels, the lifetime Tid is limited by cubic 
Dresselhaus SOI only, and as we will show later, is the same as in the completely balanced 
spin-helix case [19]. The theoretically expected values are plotted in Fig. as black lines. 



FIG. 4. Lifetime enhancement for various a//3. Lifetimes for the ID and 2D situation as 
determined from Monte-Carlo simulations for various ratios of—1.1 < a/(3 < 1.1. Data is obtained 
by fitting S z (y , t ) with a model that includes a diffusive dilution proportional to either 1/t or 
The former is used in the 2D case for \a\ ~ \/3\ (diamonds), the latter for isotropic SOI in the 2D 
case (rectangles) and for the ID case (crosses). The solid and dashed lines are theoretical curves 
(see supplementary information). The lifetime enhancement under lateral confinement is largest 
for a = 0 (arrow). For both the 2D case at \a\ = \/3\ and the ID case, the lifetime is limited by the 
same value given by the cubic SOI only. 
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The interpolation between T2d, Tcm and tid (yellow line in Fig. [3^) is in very good quantitative 
agreement with the experimental data. Although r° towards smaller w is not yet saturated, 
it is possible to project that cubic SOI will limit the lifetime. 

The lifetime enhancement achievable by channel confinement depends strongly on the 
ratio a//3. Figure [5] shows lifetimes determined by Monte-Carlo simulations for —1.1 < 
a/(3 < 1.1. For this analysis, we determined r° directly from the decay of S z (y = 0 ,t), 
see supplementary information. Without Fourier transformation one has to account for a 
diffusion factor that reduces the amplitude, in addition to an exponential decay term. The 
diffusive dilution of electrons in 2D scales with 1/t and in ID with 1 /\[t. Interestingly, the 
spins in a 2D system, however, also decay with 1 /\ft for the isotropic SOI case [H], Solid 
and dashed lines are the theoretically expected values of r° for 2D and ID spin modes (tid, 
Ad, Fpsh), as well as for the DP case (tdp). 

We find that in a narrow channel, r° does not depend on a or 3\ and is limited by cubic 
SOI (P^) only. The same limit is reached in the 2D situation at a = ft, i.e., when the system 
is tuned to the persistent spin helix symmetry. For given SOI coefficients, the maximal 
lifetime enhancement under lateral confinement in the diffusive limit occurs for the isotropic 
case (a = 0). Close to a = (3, the lifetime enhancement is small, but a reduction of diffusive 
dilution was observed [20J . 

In conclusion, we measured the evolution of a local spin excitation in a GaAs/AlGaAs 
quantum well dominated by linear Dresselhaus SOI. Because of SOI, the lateral confinement 
leads to an increased correlation between electron position and spin precession. Using a 
real-space mapping of the spin distribution, we observe a helical spin mode accompanied by 
an enhanced lifetime for decreasing channel width. The analysis in momentum space shows 
that the long-lived components decay exponentially with a minimum rate at a finite q°. Both 
the precession length and the lifetime are in quantitative agreement with theory for the 2D 
limit, the ID limit and also for the intermediate regime. The narrowest channel in our study 
still is 10 times wider than the mean-free-path (including electron-electron scattering) and 
100 times wider than the Fermi wavelength of the electrons. At those smaller length scales, 
also a reduction of the cubic SOI contribution to spin decay was predicted [21J. 

These findings illuminate an interesting path for studying spin-related phenomena. Lat¬ 
eral confinement provides a straight forward method for achieving spin lifetimes that are 
otherwise only possible by careful tuning of SOI to the persistent spin helix symmetry. This 


facilitates the use of spins in materials with stronger SOI, such as InAs or GaSb, but also in 
group-IV semiconductors, like Si and Ge. Extending the presented method to ID systems 
in the quantized limit will be relevant for the quest for Majorana fermions when combined 
with superconductors [221 - I25] . Furthermore, the results are important for transport studies 
and transistor applications [31261123 using SOI in ID or quasi-ID systems. 
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SUPPLEMENT 


Measurement Setup 

We use a time- and spatially-resolved pump-probe technique to map out spin dynamics. 
Two Thsapphire lasers are used to generate the pump and probe laser pulses at 785 nm 
and 802 nm, respectively. Pulse lengths are on the order of 1 ps and the repetition rate is 
79.1 MHz, which corresponds to 12.6 ns between two pulses. The pulses of the two lasers 
are electronically synchronized. The delay between pump and probe pulses is controlled by 
a mechanical delay line. The linearly polarized probe beam is chopped at a frequency of 
186 Hz. The polarization of the pump is modulated by a photo-elastic modulator between 
plus and minus circular polarization at a frequency of 50 kHz. The sample is located inside 
an optical cryostat at a temperature of 20 K. Both laser beams are focused onto the sample 
surface with a lens inside the cryostat. The Gaussian sigma width of the intensity profile is 
~ 1.1 /jin for both spots. The power of the pump beam is 100 /rW and that of the probe 
beam is 50 //W. The pump spot is positioned onto the sample relative to the fixed probe 
spot by a scanning mirror. After being reflected from the sample, the pump beam is blocked 
with a suitable edge filter, whereas the probe beam is sent to a detection line, where its 
polarization is monitored by a balanced photodiode bridge and lock-in amplifiers. 

Sample preparation 

A GaAs quantum well is grown on a (001) GaAs substrate by molecular beam epitaxy. 
The barrier material is Alo. 3 Gao. 7 As. Front and back Si (5-doping layers are positioned such 
that the electric field perpendicular to the quantum-well plane is very small. A sheet density 
of 3.5 x 10 15 m 2 and a transport mobility of 7.0 x 10 5 cm 2 (Vs) 1 were determined at 4 K 
after illumination by a van-der-Pauw measurement. A 5 x 5 mm 2 piece was cleaved out of 
the 2” wafer and processed with photo-lithography. Wires of variable width were etched 
by wet-chemical etching. The effective widths of the wires as given in the main text were 
determined by scanning electron microscopy images, measuring the width of the top surface. 

Theory 

The measured values of q ( y and q® in the ID limit allows the determination of a and (3. as 
described in the main text. Additionally, the knowledge of the electron density, n s , allows 
the calculation of the cubic Dresselhaus coefficient via (3 3 = —7 x n s /4 with 7 = —11 x 
10 -30 eVm 3 [2.8] . Equation ([Tj) allows the determination of D s . The structure investigated 
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is described by the following parameters. 


D s 

a 

Pi 

Pz 

0.005 m 2 /s 

-0.3 meVA 

4.9 meVA 

0.6 meVA 


The spin-dephasing time in the 2D limit for Dresselhaus fields only is given by [5| 


Dd — 


D s m 

) b 

' h A 


*2 


f3 2 + 3 0% - 2/3 2 - 


/^ 2 + P3 

8f3 2 


-1 


( 4 ) 


The spin-dephasing time in the 2D limit close to the persistent spin helix symmetry is 
given by [SI 19] 


tpsh = 2 D. 


s^r [( a _ ^) 2 + 3/3 2 ] • 


(5) 


The spin-dephasing time in the ID limit is given by da eh Eg 


Tid = 


6 ^! b 2 

ft 4 ^ 


-l 


( 6 ) 


The Dyakonov-Perel spin-dephasing time for out-of-plane spin polarization is given by 


Tdp — 


n m *2 i 1 

8 ^ (« 2 + £ 2 + /? 3 2 ) 


(7) 


The following scaling behavior is expected for the intermediate regime between the ID 
and the 2D limit 


Tim = 48t dp (q 0 w) 2 . (8) 

Real-space evaluation 

An analysis of the spin dynamics is, in principle, also possible in the real-space represen¬ 
tation S z (y, t ) and a model can be fitted to the full data. It contains the spatial variation of 
the spin mode (Bessel function in 2D, cosine in ID), a Gaussian envelope, which originates 
in the Gaussian shape of the initial spin density profile, the diffusive broadening of this 
Gaussian envelope, and an exponential decay [20]. A diffusive decay term accounts for the 
diffusive dilution of the spin density. As mentioned in the main text when discussing the 
evaluation of r° from Monte-Carlo simulations, this term is proportional to 1 j\Jt in the ID 
situation. In the 2D case, however, it is either 1/t for a ~ f3 or 1/y/t for the isotropic case, 
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FIG. 5. Fits of S z (y,t). a, Fit of a 2D real space model, as described in the text, to the data 
of the 19 — //in-channel, b, Fit of a ID real space model, as described in the text, to the data of 
the 0.7 — /im-channel. The values determined by this method for r° agree with the evaluation of 
S z (q y ,t ) in momentum-space. 

i.e., when a/3 = 0 |5j. An additional complication arises because the initial spin excitation 
has a finite spatial distribution owing to the pump-laser spot size. This can be accounted 
for by a convolution of the exact solution for a 5-peak excitation with a Gaussian function 
that is itself a convolution of the intensity profile of both the pump and the probe laser 
spots. Figure [5] shows such fits for the 19-//m and 0 . 7 -//m wires along the y direction. This 
procedure is numerically more demanding than the fits of the exponentially decaying Fourier 
components of S z (y ), and the fit parameter r° is obtained more indirectly. We therefore 
prefer to evaluate the Fourier-transformed S z . 

Fourier transformation 

To obtain the decay dynamics at specific wave-vectors, the spatial spin pattern S z (x,y) 
needs to be Fourier-transformed. If we consider a ID system defined by a channel along y, 
S z only varies along y, which allows us to apply a ID Fourier transformation: 

/ OO 

cos(q y y)S z (y,t)dy . (9) 

-OO 

In the 2D situation, the Fourier transformation in principle requires full knowledge of 
S z (x,y). Because of the almost isotropic SOI, we can assume a radially symmetric mode 
and obtain the Fourier component S z (q y ,t ) from scans of S z along y: 

POO 

S z (q y ,t)= 2nyJ 0 (q y y)S z (y,t)dy . (10) 

do 

Here, J 0 is the zeroth-order Bessel function. When the channel width w is gradually 
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FIG. 6. Comparison of radial and linear Fourier transformation Experimental values of 
t° as obtained by 2D (diamonds) and ID (circles) Fourier transformations of S z (y,t). Also shown 
are the theoretical values of t\D, toD and tim, as well as their interpolation. 

reduced, the system undergoes a transition from the 2D to the ID situation. We have 
Fourier-transformed all data sets with both methods and fitted them as described in the 
main text. Figure [6] shows r° along the //-direction determined with both methods. The 
values are very similar. In the main text, we therefore plot the 2D transformation for 
w > 15 /mi and the ID transformation for narrower channels. 

The sign of a 

While /5 can only be positive, the sign of a depends on the direction of the perpendicular 
electric field with respect to the growth direction. To determine the sign of a, we measure 
spatial spin maps also at an applied external magnetic field. Figure [7j shows such mea¬ 
surements in a channel with w = 1.7 /mi for an external magnetic field of £> ext = +1 T 
perpendicular to the wire direction. Evaluating these measurements as done in m , we can 
conclude that a < 0. Moreover, the continuous lines of constant spin orientation in these 
measurements demonstrate the helical nature of the ground mode. 

Monte-Carlo simulations 

Spin dynamics in a laterally confined 2D electron gas are calculated numerically using 
a Monte-Carlo method where the positions and spin orientations of 3 x 10 5 electrons are 
updated in time steps of 0.1 ps. Electrons are distributed on a Fermi circle and scatter 
isotropically, with the mean scattering time given by r = 2 D/vp, where vp = Tikp/m is 
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FIG. 7. Spin maps at an external magnetic field, a, S z (y,t) in a 1.7-/im-wide wire along the 
y direction at an external magnetic field of £> ex t = +1 T. The field is perpendicular to the wire 
direction. The sample temperature is 30 K. b, S z (x,t ) in a 1.7-/im-wide wire along the y direction 
at an external magnetic field of £> ext = +1 T. The field is perpendicular to the wire direction. 
The sample temperature is 10 K. The magnetic field induced additional spin precession lines with 
constant phase are tilted, showing the helical nature of the spin mode. From the opposite sign of 
the tilts in the (x,t) and the (y, £) planes, it is concluded that a is negative. 

the Fermi velocity. Each electron moves with the Fermi velocity and sees an individual 
spin-orbit field as defined in the supplementary information of Ref. m that depends on 
its velocity direction. The real-space coordinates and the corresponding spin dynamics are 
calculated semiclassically. We initialize the electrons at t = 0 all with their spins oriented 
along the 2 direction and distribute their coordinates in a Gaussian probability distribution 
with a center at x = y = 0 and a a width of 0.5 /an. Histograms of the electron density 
and the spin orientations are recorded every 5 ps, and the simulation is run until t = 5 ns is 
reached. We obtain the spin polarization at x = y = 0 versus t from the spin-density maps 
using a convolution with an assumed Gaussian probe spot size of 0.5 fi m. We determine the 
spin lifetimes r° by fitting the transients with a function proportional to 1/t x exp— t/r° 
or 1 j\ft x exp— t/r° in a window 800ps < t <4000ps, where additional spin decay is 
negligible because of the small spot sizes [19.J. For the data shown in Fig. 4, we have used 
the following parameters: D s = 0.004 m 2 /s, n s = 3.4x 10 15 cm -2 , f3i = 4.9 x 10 -13 eVm 
and /?3 = 0.6 x 10 -13 eVm. Lateral confinement was implemented by assuming specular 
scattering at the channel edges. For the ID case, w = 0.4 /rm was used. 
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